The Compound Effect of Spatial and Temporal Resolutions on the Accuracy of Urban Flood Simulation

Flood disaster is one of the critical threats to cities. With the intellectualization tendency of Industry 4.0, refined urban flood models can effectively reproduce flood inundation scenarios and support the decision-making on the response to the flood. However, the spatiotemporal variability of rainfall and the spatial heterogeneity of the surface greatly increase the uncertainties in urban flood simulations. Therefore, it is crucial to account for spatiotemporal variability of rainfall events and grids of the model as accurately as possible to avoid misleading simulation results. This study aims to investigate the effect of temporal resolutions of rainfall and spatial resolutions of the model on urban flood modeling in small urban catchments and to explore a proper combination of spatiotemporal schemes. The IFMS Urban (integrated flood modeling system, urban) is used to construct a one-dimension and two-dimension coupled urban flood model in the typical inundated area in Dongguan, China. Based on five temporal resolutions of rainfall input and four spatial resolutions, the compound effect of spatiotemporal resolutions on the accuracy of urban flood simulations is systematically analyzed, and the variation characteristics are investigated. The results show that the finer the temporal resolution is, the higher the simulation accuracy of the maximum inundated water depth. Considering the spatial resolution, as the spatial grid becomes smaller, the relative error of the maximum inundated water depth decreases, but it also shows some nonlinear characteristics. Therefore, the smaller grid does not always mean a better simulation. The spatial resolution has a greater impact on the flood simulation accuracy than the temporal resolution. The simulation performance reaches the best when the grid interval is 100 m and the rainfall input interval is 5 min, 10 min, or 15 min. Affected by other factors such as terrain slope, the simulation accuracies under different spatiotemporal resolutions present complex nonlinear characteristics. The mechanisms of the compound effect of the spatiotemporal resolutions on the model simulation and the effect of underlying surface and topography on model simulation will be the focus of in-depth exploration for the future urban flood model.


Introduction
With booming urbanization, the flood-inducing factors and hazard bearing bodies have experienced great changes in recent years. Many large and medium-sized cities around the world have suffered from frequent floods, which seriously threaten the safety of life and property of urban residents [1][2][3]. As an important basis for urban flood emergency control and risk management, the urban flood model is very important in real-time simulation, early warning, and risk assessment of floods. In recent years, under the background of Industry 4.0, the algorithms, calculation data, and computing power have greatly improved [4]. e spatiotemporal resolutions of urban flood models are becoming higher, and the decision-making is more intelligent. However, the fine-resolution simulation brings computational pressure and increases the uncertainties in the urban flood simulation. Due to the diverse underlying urban surface, the fast runoff process, and the short response time, its hydrological characteristics in urban regions show high spatiotemporal heterogeneity [5,6]. erefore, the urban flood model generally has strong resolution dependence, and the simulation accuracy is constrained by the spatiotemporal resolutions [7][8][9]. To improve the simulation accuracy, we need to conduct in-depth and systematic research on the impact of spatiotemporal resolutions on the simulation accuracy of urban flood models.
Rainfall is one of the key driving factors of urban hydrological processes and is of high spatiotemporal variability [10,11]. Variations of the rainfall spatiotemporal resolutions can affect the rainfall-runoff response time and water yield in the hydrological model [11]. erefore, it is important to simulate the hydrological response with an appropriate spatiotemporal resolution of rainfall [12,13]. A large number of studies have proved that for small urban catchments, urban flood simulation should use the rainfall data on at least 1-15 min time resolution and 100-1000 m spatial resolution [14][15][16]. Bruni et al. [17] analyzed the relationship between the urban catchment area and the required spatiotemporal resolutions of rainfall data and showed that for small urban catchments with less than 700 hectares, the urban hydrological simulation requires the rainfall data with at least a 5 min time resolution and a 1.7 km spatial resolution. In addition, the sensitivity of different hydrological models to the spatiotemporal resolutions of rainfall differs significantly, with the physical model being more sensitive than the conceptual model [18][19][20]. Aronica et al. [21] found that the Storm Water Management Model (SWMM) is more sensitive to rainfall temporal resolution than hydrological parameters. Meselhe et al. [22] compared the HMS (hydrologic modeling system) conceptual model with the physics-based hydrological model MIKE SHE. ey found that the latter is more sensitive to rainfall temporal resolution. Gires et al. [23] found that the 1D/2D (one-dimension/two-dimension) coupled model, Multi-Hydro, is more sensitive to rainfall variability than the simpler 1D model. Many other studies also found that in small urban catchments, the time-scale variation of rainfall data has a greater impact on urban hydrodynamic models [24,25]. To sum up, for a densely built and highly impermeable urban catchment, the output of the urban flood model is very sensitive to the rainfall spatiotemporal resolutions, and the bias of output increases obviously as the resolutions decrease.
Most urban flood models are distributed hydrological models based on grid data. e hydrological process is also sensitive to the spatial distribution of the underlying surface of the watershed [9]. e computational efficiency of the model and the accuracy of the simulations are often affected by the spatial resolution of the grid and the accuracy of input data [26,27]. In the early 1960s, the importance of the spatial resolution of input data is recognized by scholars [28,29]. Since then, many scholars have studied the spatial resolution of hydrological models. ey found that the conclusions are different in different study areas. e model of high spatial resolution can lead to the systematic underestimation of peak flow [30,31]. Ichiba et al. [27] found that increasing spatial resolution can reduce peak flow and total flow. e effect of spatial resolution of the grid on model performance is nonlinear, and higher mesh accuracy does not necessarily lead to better simulation. Zhang and Montgomery [32] found that a spatial resolution of 10 m can greatly improve the simulations than 30 m and 90 m, while that of 2 m or 4 m only improves the model slightly. e above studies show that the difference in spatial resolutions will change the loss of underlying surface information and the complexity of the surface runoff between adjacent grids, which has a nonlinear impact on the simulations. erefore, in the modeling process of the urban flood, the appropriate spatial resolution should be selected according to the comprehensive analysis of the calculation characteristics of the model, the characteristics of the underlying surface of the study area, and the accuracy of the input data.
To sum up, there have been extensive studies on evaluating the accuracy of urban flood simulation unilaterally from the temporal resolution of rainfall data or spatial resolution of the grids. ese studies show that the effects of temporal resolution and spatial resolution on the simulation accuracy are obvious and nonlinear. So, the compound effects of grid spatial resolution and rainfall temporal resolution on the simulation accuracy should be more complicated. However, there are few studies on this aspect. To this end, we take the typical flooded area in Dongguan, China, as the study area (Section 2.1) and collect topographic and sewer network data (Section 2.2). e distributed urban flood model IFMS Urban (Section 3.1) is used to construct the coupling model of the 1D urban drainage network model and 2D surface hydrodynamic model in the study area, and the model is validated by using the high temporal resolution rainfall data and observed historical inundation events (Sections 2.2 and 3.3). en, we explore the effects of rainfall temporal resolution and model spatial resolution on the maximum inundation water depth and the submerged water depth hydrograph based on five temporal resolutions of rainfall input and four spatial resolutions(Sections 4.1 and 4.2). Furthermore, we analyze the compound effects of different spatiotemporal resolutions on the accuracy of urban flood simulations (Section 4.3). Conclusions can be found in Section 5.

Study Domain.
e study site, Guancheng District of Dongguan City, is located on the south-central east coast of the Pearl River estuary, Guangdong Province, China (within E119°31′-114°15′, N22°39′-23°09′). It covers approximately 13.31 km 2 (Figure 1(a)). e study area represents a typical urban area with a density of human structures such as houses, commercial buildings, and roads. When the city is hit by heavy rainfall, the topographic characteristics of the study area make the flood converge to the middle from the sides south and north and finally discharge into the Dongyin Canal in the east.

Topographic Data.
e basic geographic data provided by the Urban Planning Bureau of Dongguan City (UPBDC) included a remote-sensing image and a digital elevation model (DEM). e resolution of the remote-sensing and 2 Computational Intelligence and Neuroscience DEM are 5 m and 0.5 m, respectively. e former was used to distinguish the land use types of the underlying surface, and the latter was used to calculate the elevation and slope of the 2D grid.

Sewer Network Data.
e sewer network data (Figures 1(b) and 1(c)) that were obtained from UPBDC is mainly based on the combined system. Most of the pipelines do not meet the design return period of one year, and most of the pipelines are less than 1000 mm in diameter, which is easy to cause waterlogging in the lower terrain.

Rainfall Data.
In this research, two rainstorm events that occurred on August 30, 2018, and May 7, 2015, were used for model calibration and impact study of temporal and spatial scale, respectively. e maximum 24-hour cumulative rainfall of the two rainstorm events exceeded 50 mm, resulting in waterlogging in the study area. ree rainfall stations, which collect rainfall data at an interval of every 5 min, are located in and around the study area (Figure 1(c)). e rainfall of 10 min, 15 min, 30 min, and 1 h temporal scale was accumulated from the rainstorm of May 7, 2015.
Meteorological Bureau of Dongguan City provided the corresponding rainfall records.

Observed Historical Inundation Events.
Generally, it is difficult to calibrate and validate the urban inundation model due to a lack of detailed observation of inundation events [2]. erefore, in many cases, either a partial calibration/validation of the model or an indirect validation/ verification based on testimonial reports is sought [33,34]. In this research, the process of inundation depth in the inundation area in Yonghuating, Dongcheng Road West, and Dongzong Road on August 30, 2018, was observed for the calibration and verification of the urban flood model. e data on maximum inundation depth and distribution of inundation area on May 7, 2015, were collected and released by the Water Authority of Dongguan City.

Methodology
3.1. IFMS Urban. IFMS Urban couples the SWMM with a 2D surface hydrodynamic model that is conducted based on an adaptive grid and finite volume method, which can automatically identify the region with a large parameter Computational Intelligence and Neuroscience gradient and the boundary between dry and wet, adjust the grid size, accurately simulate the dynamic change of water flow propagation, and can be applied to the actual flood simulation [35]. e IFMS Urban efficiently calculates the urban flood process in the complex urban area with frequent waterlogging and achieves good calculation accuracy, which provides technical support for this paper.

1D Urban Drainage Network
Model. e urban drainage system consists of water inlets, drainage pipes, drainage pumping stations, and river channels at the outlets of the pipe network. e 1D model can use three methods of dynamic wave method, kinematic wave method, and steady flow method to calculate the drainage pipeline confluence. e governing equation is specifically the following formulas: where Q is the discharge, A is the discharge section area, t represents the time, and x represents the distance.
where H represents the water depth, g represents the gravitational acceleration, and S f is the slope gradient of friction resistance.

2D Surface Hydrodynamic Model.
In order to establish a special well-balanced scheme technique for dealing with source term due to bottom topography constructed, this paper develops a well-balanced Godunov-type scheme of the second-order accuracy for 2D shallow water equation with mesh. As above, the MUSCL method is used to reconstruct the variable values on both sides of the unit interface U L/R i+1/2 and the Roe format is selected to solve interface flux in the evolutionary step. Regarding the discretization of the source term, the bed slope term is discretized by characteristic classification, and the resistance source term is discretized implicitly. e 2D shallow water equation of depth-averaged can be abbreviated as follows: where h represents the water depth, u is the flow velocity of x-direction, v is the flow velocity of y-direction; S x and S y are the source terms.

Algorithm of Coupled
Model. e 1D urban drainage network model and 2D surface hydrodynamic model are coupled by calculating the exchange water volume, which is substituted into their respective model for calculation and update to the next step. e exchange water volume can be calculated through the following equation: where H Surface is the head of land surface, H node is the head of drainage pipeline, M is the discharge coefficient, W crest is the width or perimeter of manhole, and H g is the surface elevation.

Evaluation of Modeling
Results. e calibration process has been evaluated using error indicators including Nash-Sutcliffe efficiency coefficient (NSE), maximum inundation depth relative error (RE P ) , and maximum inundation depth appearance time absolute error (AE T ) for the simulated and observed values. e error indicators are formulated as in the following equations: where q obs i is the observed value of an event i, q sim i is the simulated value of an event i, N is the number of observed values, and q obs is the average of observed values.
where q obs P is the observed value of maximum inundation depth, q sim P is the simulated value of flood peak maximum inundation depth.
where T obs P is the observed value of maximum inundation depth appearance time, and T sim P is the simulated value of maximum inundation depth appearance time.
In order to compare the results from different rainfall temporal resolution and 2D grid scale and also to compare the effect of the composition of the rainfall temporal resolution and 2D grid scale, different measures were used. In addition to the common error indicators such as RE P and AE T , we also used the coefficient of determination (R 2 ) to compare the correlation of two inundation depth series. e closer R 2 to 1, the higher the correlation between the two inundation depth series. R 2 is specifically as follows: where Y A is the value of the inundation depth series A, Y A is the mean value of the inundation depth series A, and Y B is the value of the inundation depth series B.

Model Validation
e model is set up based on the drainage network. e study area is generalized into 999 pipes, 1011 manholes, and 12 outlets. e parameters for the model include the section size of the drainage network, the two-dimensional (2D) grid elevation, the subcatchment slope, and the 2D grid pervious surface ratio. e section parameters of the drainage network are obtained from the actual survey data. e 2D grid elevation and the subcatchment slope are determined based on the DEM data. e 2D grid pervious surface ratio is extracted from the remote sensing images over the study area. e empirical parameters of the model include the impervious area Manning's roughness, the depression storage for the pervious area, and the infiltration related parameters. e empirical parameters are calibrated according to the SWMM user manual and the hydrogeological characteristics of the study area, as shown in Table 1. e rainwater in the study area flows through the pipeline by gravity and is finally discharged into the Dongyin Canal. According to historical data, the water level of the Dongying Canal is low and does not affect the outflow of the pipeline. erefore, the model boundary outflow condition is set to free outflow. e modeling time step is set to 20 s, and the total simulation time is 24 h. e setting can ensure that the pipeline network system has no inflow. e accumulated water in the pipeline has been emptied when the simulation ends under each input rainfall condition.

Model Validation.
In this section, we validate the IFMS Urban model by using the rain and flood event on August 30, 2018, as an example. e rainfall poured heavily in the study area from 11:00 to 16:00, within which rainfall amount accounted for over 80% of the total in the day. us, a 5 h storm event spanning from 11:00 to 16:00 was used as model rainfall data.
A summary of the simulations performed is given in Figure 2 and Table 2. e model can well predict the flood occurring and receding at the typical waterlogging points of Yonghuating, Dongcheng West Road, and Dongzong Road. ese results confirm that the simulation of flood processes is reasonably well, with NSE for all waterlogging points exceeding 0.75. e RE P and AE T ranged approximately from 11.11% to 17.50% and from 5 min to 12 min, respectively. All of them are within the allowable error range required by the Standard for Hydrological Information and Hydrological Forecasting of China. Furthermore, the results of Dongcheng West Road are better than other waterlogging points, with NSE, RE P and AE T being 0.86, 11.11%, and 5 min, respectively. In conclusion, the above results show that the parameters of the model are set reasonably, and the model has good applicability in the study area. e model is reliable and is suitable for the subsequent analysis. Table 3 and Figure 3 show the maximum inundation water depth and RE P at five typical waterlogging points under five temporal resolutions. With the increase of rainfall temporal resolution, the maximum inundation water depth gradually decreases and RE P increases. However, the effect of the rainfall temporal resolution on the maximum inundation water depth is small on the whole. Increasing the time step of rainfall input from 5 min to 30 min, the variation range of the maximum inundation water depth at the five waterlogging points is within 5 cm, while the variation range of RE P is less than 8%. RE P at waterlogging points based on the 1 h temporal resolution is larger than the other four temporal resolutions, with that in Dongcheng West Road being the largest (13.73%). e errors of the maximum inundation    Computational Intelligence and Neuroscience water depth at the waterlogging points increase with the temporal interval of rainfall input, but the overall increase is not large. e occurrence time of the maximum inundation water depth (referred to as the peak time of inundation water depth) at waterlogging points has a high correlation with the temporal resolution of rainfall input. As shown in Table 4, with the increase of the temporal resolution of rainfall input, the time between the peak time of inundation water depth at waterlogging points and the rainfall peak time (referred to as the lag time) also increases. e lag time of the 10 min and 15 min rainfall input intervals at the five waterlogging points is short (only 0-3 min). e lag time of the 30 min and 1 h rainfall input intervals changes a lot. Compared with the 5 min time interval, the lag time is, respectively, extended by 7-13 min and 23-36 min.

Effect of Rainfall Temporal Resolution.
With the increase of rainfall input interval, the correlation between the submerged water depth and the rainfall at waterlogging points decreases. e deformation degree of the submerged water depth hydrograph at waterlogging points becomes larger. Taking the waterlogging point, Yonghuating, as an example (as shown in Figure 4), the rainfall peaks twice at about 1 h and 2 h. Correspondingly, there are two obvious waterlogging processes in 1-2 h and 3-4 h under the 5 min, 10 min, 15 min, and 30 min rainfall temporal resolutions at Yonghuating, and the results have a good correlation with the rainfall observation. e submerged water depth hydrograph under the 1 h rainfall temporal resolution has only one obvious waterlogging process. Compared with the waterlogging processes under other temporal resolutions, the maximum inundation depth is lower and the corresponding occurrence time is later. Using the coefficient of determination (R 2 ), this study further analyzes how well the shape of submerged water depth hydrographs under 10 min, 15 min, 30 min, and 1 h rainfall input resolutions at five typical waterlogging points match that of the hydrograph under the 5 min resolution ( Figure 5). With the increase of temporal interval, the coefficient of determination decreases, and the difference between submerged water depth hydrographs increases. e      Computational Intelligence and Neuroscience shape for the 10 min and 15 min resolutions is in good agreement with the submerged water depth process under 5 min resolution, and R 2 are all higher than or equal to 0.97. However, there are certain differences of hydrographs between 1 h and 5 min resolutions, and the R 2 is only 0.69-0.81.

Effect of the 2D Spatial Resolution of Model Grids.
e study area is used as the meshing area to generate four kinds of mixed 2D meshes of triangles and quadrilaterals with spatial scales of 50 m × 50 m, 100 m × 100 m, 150 m × 150 m and 200 m × 200 m. According to the remote sensing images, we analyze the land use type of the underlying surface and set the impervious area and roughness of the simulations under four spatial resolutions. e grid elevations for the four spatial resolutions are all obtained by inverse distance weighted interpolation based on a DEM with a resolution of 5 m.
As shown in Table 5 and Figure 6, the effect of grid spatial resolution on the model simulation is nonlinear. It means that the smaller spatial resolution does not correspond to the higher accuracy of model outputs. With the increase of grid spatial interval, the maximum inundation depth decreases and the absolute error (AE) of the maximum inundation depth increases on the whole. When the grid spatial resolution is increased from 50 m to 150 m, the maximum inundation depth decreases, and its variation range is small. All AE is within 4 cm and RE P is less than or equal to 6%, while those of the 200 m grid are larger than the other resolution simulations. e maximum inundation depths with a spatial resolution of 50 m at the five typical waterlogging points are all overestimated. Except for Yinshan Street, RE P of the 50 m grid is higher than that of the 100 m grid at other four waterlogging points.
ere are significant correlations between the slope and the grid spatial resolution of the catchments of waterlogging points and the hydrological characteristics. As the spatial resolution becomes larger, the calculation accuracy of the waterlogging point with large slope of catchment area decreases obviously. e DEM and flow direction of main roads in the study area and around each waterlogging point are shown in Figure 7. On one hand, the three areas of Yinshan Street, Dongzong Road, and People's Park are located in the low-lying central of the study area, where the rain and floods in the eastern, northern, and southern parts of the study area are concentrated. But due to the flat terrain from the Yinshan Street to the People's Park, the converged rain and floodwater cannot be discharged in time, resulting in waterlogging in the three areas. When the grid resolution is increased from 50 m to 200 m, the simulation results at Yinshan Street, Dongzong Road, and People's Park are less affected by the grid size. All AE is less than 3 cm and RE P is less than 13%. On the other hand, small-scale rain and floods are gathering in the mountainous southeast of the study area, mainly located in Yonghuating and Dongcheng West Road. e slope in the catchment areas is larger than that in other waterlogging areas. When the grid resolution is increased from 50 m to 150 m, AE and RE P in Yonghuating and Dongcheng West Road have a small variation range, while AE and RE P at the 200 m grid have increased obviously, reaching 7.07-27.26 cm and 23.6-45.43%, respectively. In the urban flood model, the spatial grid is the basic unit of the runoff calculation of the model. e slope and other properties of the underlying surface have a direct impact on the runoff calculation. In this study, most of the study area is continuous urban hard ground or green space. With the increase of the spatial grid interval, the generalization of the terrain and underlying surface information is greater, and variation details of the urban terrain and Computational Intelligence and Neuroscience underlying surface are covered up. erefore, the simulations at Yonghuating and Dongcheng West Road where there is larger slope are more affected by the spatial resolution. As shown in Figure 8, the effect of spatial resolution at typical waterlogging points is greater than that of rainfall temporal resolution. For the grids with different resolutions, with the increase of rainfall temporal resolution from 5 min to 30 min, the maximum inundation depth at typical waterlogging points has a small variation range, and RE P is less than 3%. e effects of the 5 min, 10 min, and 15 min temporal resolutions are very close. For different rainfall temporal resolutions, with the increase of grid spatial interval, RE P at typical waterlogging points greatly changes with nonlinear characteristics. For the waterlogging points of Dongcheng West Road and Yonghuating with a large slope, large grids have a larger influence on the simulation results, while for the flat Dongzong Road and People's Park, small grids have a larger influence on the simulation results.

e Compound Effect of Rainfall Temporal Resolution and
In urban flood simulation, the finer spatiotemporal resolution does not mean more accurate simulation results. For the 50 m grids, the maximum inundation depth is overestimated under different rainfall temporal resolutions. e finer the temporal resolution is, the greater the overestimation is. At the five typical waterlogging points, RE P with the spatiotemporal resolutions of 5 min and 50 m is relatively large, ranging from 1.8% to 10%. In contrast, the water depth range that decreases with the increase of rainfall temporal resolution is small. It is not enough to compensate for the overestimated water depth range under the highprecision grid and RE P is still larger than that of the 100 m grid. For the grid of 100 m, the overall simulations are better than those of 50 m and 150 m. e 5 min × 100 m, 10 min × 100 m, and 15 min × 100 m resolutions can achieve     Computational Intelligence and Neuroscience the best simulations at different typical waterlogging points, and their RE P ranges from 0.3% to 6.1%. For the grid of 150 m, the maximum inundation depth of typical waterlogging points is underestimated. With the increase of rainfall temporal resolution, the maximum inundation depth decreases and the RE P increases. By judging the compound effect of rainfall temporal resolution and grid spatial resolution on the simulations, we find the simulation results are better with the 5 min, 10 min, or 15 min rainfall temporal resolution and the 100 m grid spatial resolution.

Conclusions
is study utilized the IFMS Urban model that couples SWMM and 2D hydrodynamic model to analyze the response of the flood process of urban to different spatiotemporal resolutions. e typical waterlogged area, Guancheng District of Dongguan City, China, was selected as the study area. e variation characteristics of maximum inundation water depth and inundation process were investigated based on five rainfall temporal resolutions and four grid spatial resolutions. e main conclusions are as follows: With the increase of rainfall temporal resolution, the accuracy of the simulated maximum inundation depth decreases, the time interval between the occurrence of the maximum inundation depth and the peak time of rainfall is prolonged, and the shape of the submerged water depth hydrograph changes greatly, which means the inundation process correlation with the rainfall process is decreased. Compared with the temporal resolution of 5 min, 10 min, 15 min, or 30 min, the 1 h temporal resolution has a greater influence on urban flood simulation. For the effect of grid spatial resolution, the effect of grid resolution on urban flood simulation is nonlinear. With the increase of grid interval, the maximum inundation depth decreases and RE P presents an overall increasing trend. For the fine 2D grids, the maximum inundation depth is overestimated, and the simulations at the waterlogging points with larger slopes are more affected by the spatial resolution.
In this study, the grid spatial resolution has a greater effect on the simulations at typical waterlogging points than the rainfall temporal resolution. It is not the case that the finer spatial and temporal resolutions are, the higher the accuracy of the simulation will be. e simulations perform better with the 5 min, 10 min, or 15 min rainfall temporal resolution and 100 m grid spatial resolution. For the area with the large slope in the catchment areas, using a 2D grid with a smaller grid spatial interval can obviously improve the computational accuracy of the model.
In this study, we have not considered the influence of terrain slope and vertical structure on the urban flood model. In the future, we will focus on the mechanisms of the compound effect of the spatiotemporal resolutions, and on how the model accuracy is affected by the type and slope of underlying surface.

Data Availability
e rainfall records are provided by the Meteorological Bureau of Dongguan City. e data of maximum inundation depth and distribution of inundation area on May 7, 2015, were collected and released by the Water Authority of Dongguan City. e basic geographic data provided by the Urban Planning Bureau of Dongguan City included a remote-sensing image and a digital elevation model. e sewer network data was provided by the Urban Planning Bureau of Dongguan City.

Conflicts of Interest
e authors declare that they have no conflicts of interest to report regarding the present study.